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• and Marcus 0. Creager 


SUMMARY 


The equilibrium temperature and heat-transfer coefficients for 
transverse cylinders in a high-speed stream of rarefied gas were 
measured over a range of Knudsen numbers (ratio of molecular -mean-free 
path to cylinder diameter) from 0.025 to 11.8 and for Mach numbers from 
2.0 to 3*3* The range of free-stream Reynolds numbers was from 0.28 to 
203. 


The models tested were 0.0010-, 0.0050-, O.O 3 O-, O.O 5 I-, O.O 8 O-, 
and 0 . 126 -inch-diameter cylinders held normal to the stream. 


For the case of high-speed flow of a rarefied gas about cylinders, 
the data indicate that fully developed free-molecule flow first occiirs 
at Knudsen numbers of approximately 2.0 and the temperature -recovery 
factor depends primarily upon the Khudsen number. For Knudsen numbers 
greater than 0.2, the teraperatiire— recovery factor exceeds a value of 
unity even though free-molecule flow may not be fully developed. 


/ 


Over the range of conditions covered in the present tests, the 
Nusselt number is a function only of the Reynolds n\miber if the vis- 
cosity and thermal conductivity are based upon stagnation temperature 
and the density based upon free-stream conditions. For free-molecule 
flow, the heat -transfer data are well correlated by the theory with the 
accommodation coefficient equal to O. 9 O. 


INTRODUCTION 


The subject of heat transfer to or from bodies in a high-speed 
rarefied gas stream has become of importance in aeronautics . The recent 
practical realization of missile flight at very high altitudes has 
stimulated interest in the subject. Also, in certain aspects of 
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wind-tunnel testing^ notably turbulence studies in supersonic wind 
tunnels by means of hot-wire anemometers^ the effects of slip may not 
be negligible. (See^ e.g.^ reference 1.) 

The various phenomena encountered in high-speed flow of gases are 
often described in terms of the Reynolds number and the Mach number of 
the flow. The former parameter may be considered to be a measure of 
the effect of viscosity and the latter a measure of the effect of 
compressibility on the flow field. 

When considering the flow of rarefied gases ^ however^ an additional 
parameter^ the Knudsen number^ becomes important. The Knudsen number 
is defined as the ratio of the mean-free molecular path to a character- 
istic body dimension. This parameter^ which is a measure of the degree 
of gas rarefaction, may also be construed as defining the importance of 
the effect of molecular motions on the flow field. When the Knudsen 
number is small (less than, say, O.OOl), the effect of molecular motions 
on the flow is negligible and, in this regime, the gas can be treated 
as a continuous medium. On the other hand, for large values of the 
Knudsen number (of the order of 10), the effect of molecular motions is 
all-inportant and the phenomena which occur can be completely described 
in terms of the motions of individual molecules - in a statistical 
sense, of course. By virtue of the well-established laws of molecular 
motions, this so-called "free -molecule" regime is readily susceptible 
to analysis and has been thoroughly examined (references 2 through 8) 
with respect to both aerodynamic and thermodynamic phenomena. Experi- 
mental work in this regime with respect to external flows (as opposed 
to flow in tubes) is, however, quite meager. The pioneer work of 
Epstein in 1924 (reference 3) '^^as concerned with the measurement of the 
drag of spheres; the tests were confined to very low-speed motions. 
Recently, the drag and equilibrium tenperature of a cylinder in high- 
speed free -molecule flow were reported in reference 8. 

The regime of intermediate Knudsen numbers (O.OOl < K < lO) is 
difficult to deal with analytically and, at this time, little success 
has been attained in attempts to analyze this portion of the field of 
fluid mechanics. Some progress has been made, however, by Chapman and 
Cowling (reference 9) > Schamberg (reference lO) , and Grad (reference ll) . 

The two general types of approach to the problem used by the above 
authors have been described by Uhlenbeck (reference 12) . Briefly, the 
first approach consists of applying correction terms to the Navier- 
Stokes continuum equations and boundary conditions, the first approxi- 
mation resulting in a set of higher-order equations - the Burnett 
equations. The other approach is to attempt to solve the Boltzmann 
equation for the molecular velocity distribution, which, in effect, 
results in obtaining correction teims to the Maxwell velocity distri- 
bution equation. Both of these methods are essentially perturbation 
methods . In the first approach solutions are obtained in terms of 
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power series of the Knudsen number^ while in the second method solutions 
are obtained in terms of the reciprocal of the Knudsen number. Both 
methods will fail for Knudsen numbers of the order of unity. 

As in the case of free -molecule flows, little experimental work 
has been done in the intermediate regime for the case of high-speed 
external flows. However, Kane (reference 13) has reported results of 
drag tests on spheres in this regime. 

The general puip)ose of the research described in this paper was to 
study the heat-transfer processes that occur between a body and a high- 
speed rarefied gas stream. More particularly, it was hoped that the 
experimental data would yield two important results: first, the proper 

grouping of variables to obtain correlation of heat-transfer data in 
the analytically intractable "slip -flow” regime intermediate between 
the continuum and free-molecule regimes; and second, to define approxi- 
mate values of the parameters which delineate the three regimes for the 
particular configuration used in these tests. 

Experimental data were obtained on the heat transfer from circular 
cylinders in the three regimes of free-molecule, slip, and continuum 
flows . The data are compared at small values of Knudsen number 
with the data of reference 1 and, for large values of Knudsen number, 
with the free-molecule analysis of reference 8. 


NOTATION 


Cp specific heat at constant pressure, foot-pounds per slug, ^ 

d cylinder diameter, feet 

f(s) dimensionless function of s defined by 
f(s) = Zi(s2 +3) + 22^32 + 1^ 

g(s) dimensionless function of s defined by g(s) = 3(Zj^ + Zg) 

g gravitational acceleration, 32.2 feet per second squared 

h heat-transfer coefficient, foot-pourids per square foot, 

second, 

Iq modified Bessel function of first kind and zero order 

Ii modified Bessel function of first kind and first order 

K Knudsen number dimensionless 

k Boltzmann constant, 5-66 X 10”^'* foot-pounds per °F per molecvile 
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I 

M 

Nu 

P 

Pr 

Q 

r 

R 

Re 

s 

T 

Tc 

Te 

U 

V 


a 


mean-free molecular path, feet 
Mach number f ionless 


Nusselt number [ ~ ) , dimensionless 

\ 

static pressure, pounds per square foot 
Prandtl number dimensionless 

V 

heat transferred from body, foot-pounds per square foot, second 

temperature-recovery factor, defined by Te = T ^1 + r M‘ 
dimensionless 

gas constant, foot-pounds per pound, ^ absolute 
Reynolds number f dimensionless 

molecular speed ratio ( ^ dimensionless 

\ 

free-stream static temperature, ^ absolute 
temperature of the cylinder, ^ absolute 
equilibrium temperature, ^ absolute 
stream mass velocity, feet per second 
acoustic speed ( ^7gRT) , feet per second 


most probable molecular speed (v^ 2gRT) , feet per second 

^8gRT 


mean molecular speed 


), feet per second 


dimensionless function of s defined by 


dimensionless f-unction of s defined by 


exD 


= Jts2 


Io( ^ ) + I- 


1 y 

exp ( - 


2 y • 2 

accommodation coefficient, dimensionless 
ratio of specific heats, dimensionless 
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K 

X) 

p 


thermal conductivity, foot-pounds per second, square foot, 
^ per foot 


viscosity, pound-seconds per square foot 


kinematic viscosity 


(0 


, square feet per second 


gas density, slugs per cubic foot 


Subscript 


10 refers to the fact that the viscosity and thermal conductivity 
were evaliiated at stagnation temperature while the density 
was evaluated at free-stream conditions 


AJJTALYSIS 

Some General Relations Between Variables 


Before proceeding with the analysis, it is desirable to set down a 
few relations which are used throughout. The Knudsen number, K = Z /d, 
can be related to the Mach and Reynolds numbers as follows: 


iL = 

Re 


( 1 ) 


From reference page 1^7, the viscosity of a rarefied gas composed of 
hard elastic spheres and having a Maxwellian velocity distribution can 
be expressed in terms of the density, mean molecular velocity, and mean- 
free molecular path, = 0A99P v I . Thus, very closely. 


iL 

Re 



K 


( 2 ) 


The molecular speed ratio s, which is a primary variable in the 
analysis of free-molecule flows, is defined as s = U/vm. Since 
vm = 2gRT and Va = </t^T, we have a relation between s and M 



s 


M 


( 3 ) 
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Heat Transfer From Cylinders in Free -Molecule Flow 


Heat transfer from a circular cylinder oriented transversely to a 
gas stream has been treated, for the case of free -molecule flow, in 
reference 8. It is useful, however, to examine the final equations 
derived therein with a view to putting them in such form as will allow 
plotting of the experimental data on a basis comparable with continuum 
data • 


The eqiiation for heat transfer from a cylinder to a diatomic gas 
stream is shown in reference 8 to be (neglecting radiation terms) 


^ g(s) - f(s) 


pv^a 


Q = 0 




One can define a heat-transfer coefficient per unit area as 



(5) 


where Tg is the equilibriiom teiiperature ass\jmed by the cylinder in 
the absence of heat transfer. The quantity Te can be obtained from 
equation (4) by setting the term Q equal to zero in which case Tg 
equals Tg- Thus, 


Te 



( 6 ) 


Then, substitution of eqmtions (4) and (6) into equation (5) yields 


pvma g(s) 


( 7 ) 


If use is made of the relation p/T = gpR, equation (7) can be rewritten 
as 


Hu = -S-fJ— Re Pr (8) 

Substituting numerical values for the constants, for the case of 
nitrogen gas, we have, 

Nu = 0.0264 a Re Pr 


s 


(9) 
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The term g(s)/s is nearly constant for values of s greater than 2 
and approaches a value of 6»/~n asymptotically. The limiting values 
are 


g(s) 

s 


s-> o 


s 


g(s) 


S—> 00 


6 ^/n 


Then for values of s > 2 we have, very nearly, 

Nu = 0.283 a Re Pr 


( 10 ) 


For gases, the Prandtl number, Pr, is nearly constant over a wide 
temperature range and, for diatomic gases, has a value of about 0.72 
at 0^ F . Substitution of this value for Pr into equation (lO) yields 
the final equation for heat transfer from a transverse cylinder in a 
high speed (s > 2) free -molecule flow field 

Nu = 0.204 a Re (ll) 


DESCRIPTION OF EQUIPMENT AND EXPERIMENTAL TECHNIQUE 

The Wind Tunnel 


The Ames Laboratory low-density wind tunnel, in which the present 
tests were conducted, is an open-jet, nonreturn -type tunnel and is 
described in detail in reference 8. A sketch of the major components 
of the wind tunnel is shown in figure 1 . 

The nozzles used in the present test, however, differed from those 
described in reference 8 and for that reason will be described in detail 
here. The nozzles used were axially symmetric, l-l/2 inches in exit 
diameter, and had porous walls to allow part of the boundary layer to 
be removed by suction. The nozzle contours were determined by calcu- 
lating the shape of the inner core (assumed inviscid) by the method of 
characteristics, as outlined in reference 15; and applying a boundary- 
layer correction to the contour. The boundary -layer correction was 
calculated using the von Karman momentum- integral theory given by Schaaf 
in reference l6. The correction actually applied to the nozzle contours 
corresponded to that calculated for the highest pressure level used. 

The nozzles were constructed by stacking 2-inch-diameter, thin, 
metal shims to the desired length. The stack of shims was held together 
by four through-bolts and was machined to the calculated internal 
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contour. The nozzle walls were made porous by discarding every other 
solid shim and replacing it with small spacers around each of the four 
through -bolt s . The shim thickness was varied along the length of the 
nozzle in order to obtain a reasonably uniform suction velocity normal 
to the wall. The nozzle boundary -layer removal was accomplished by 
connecting a chamber around the nozzle to one of the main vacuum pumps. 
The porous nozzle is shown diagrammatically in figure 2. 

The use of boundary -layer removal resulted in a nozzle which pro- 
duced an essentially constant Mach number over the center third of its 
diameter in contrast with the solid nozzles described in reference 8 
which produced essentially a parabolic Mach number distribution. Also 
the use of boundary -layer removal permitted the pressure level of the 
wind tunnel to be changed by a factor of 2-l/2^ while the Mach number 
varied only by about ±2 percent . A typical Mach number distribution 
obtained with and without boundary -layer suction is shown in figure 3. 

The three nozzles used in i:he tests produced Mach numbers of 2.00 
11-1/2 percent^ 2.50 ±2 percent^ and 3«15 12 percent over a range of 
pressure levels from 80 to 200 microns of mercury absolute in the center 
half -inch of the jet. 


Models Tested 


The models tested were right circular cylinders held normal to the 
air stream by copper rods fastened to the ends of the cylinders. Six 
models were used in the tests; their diameters^ test lengths^ and types 
of construction are given in the following table: 


Model 

no. 

Material 

Diameter 

(in.) 

Test 

length 

(in.) 

Remarks 

1 

Pt-Ni alloy 

0.0010 

0.434 

Solid 

2 

Pt 

.0050 

.480 

Solid 

3 

Fe - Constantan 

.030 

.50 

Butt -welded 





thermocox:^le 

k 

Pt film on 





Pyrex glass 

o 

CO 

O 

.530 


Solid 

5 

Ni 

.051 

.380 

j 

Ni wire woiind 






1 on ceramic 

6 

m 

.126 

.368 


i tube 
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Models 1 through k had guard heaters wound on the copper support 
rods to allow the temperature of the ends of the models to be adjusted 
in order to prevent end losses. Thermocouples were soldered to each 
end support in order that the model end temperature could be measured. 
Figure k shows the arrangement of the model and its supports in the 
stream. Models 5 and 6 were constructed by wrapping 0.003-inch-diameter 
insulated nickel wire on a ceramic tube. Voltage leads were soldered to 
each end of the test length and brought out through the ceramic tube to 
the ends of the model. The model was then dipped in a thin solution of 
enamel to obtain a smooth surface. These models were made long enough 
to completely span the nozzle. 

Model 3 ^as a butt-welded iron-constantan thermocouple with the 
junction located in the center of the test length. 

Model k was constructed by coating a Pyrex glass rod with a very 
thin layer of platinum. The model was moimted as shown in figure k by 
soldering the copper rods directly to the platinum coating. 

The temperature of all the models^ except number 3^ was determined 
from their resistance which was measured by comparing the voltage drop 
across the model to the voltage drop across a standard resistor connected 
in series with the model. These voltages were measured with a direct- 
current potentiometer and galvanometer having a least count of 1 micro- 
volt. The temperature-resistance characteristics of each model were 
determined by calibrating the model in a constant tenperature bath over 
the range of temperatures encountered in the tests. The temperature of 
model 3 was determined by measuring, with a potentiometer, the thermo- 
couple voltage generated at the junction. 


Test Procedure 


An impact -pres sure tube, a static-pressure tube, and the test 
model were mounted on a common support which could be moved through the 
test area of the jet. The wind tunnel was then started and the impact 
and static pressures were measured over the test area of the nozzle. 

The static-pressure tube has been described in reference 8 and the impact- 
pressure tube was identical with the 0.150-inch-diameter impact-pressure 
tube described in reference 17 . Pressure-rmeasuring techniques are also 
described in reference 8. The impact-pressure corrections due to viscous 
effects, given in reference 17 , were used in order to calculate the 
stream Mach number. 

After the stream conditions had been determined, the test model was 
lowered into the stream and the equilibrium temperatures of the model 
test length and ends were measured. In the case of models 1 and 2, the 
guard heaters were used to raise the model end temperatures until they 
were equal to the indicated test-length temperature. Heat was not 
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supplied to the ends of the other models for the measurement of equi- 
librium temperature as, for the most part, the equilibrium temperatures 
were, on the average, IQO F below the stagnation temperature of the 
stream. 

After the equilibrium ten^'erature of the model was measured, the 
model was heated by passing a direct current from a battery through the 
model. For models 1, 2, and the end temperatiires were adjusted until 
there was no temperature difference between the ends and the model test 
length. When this point was reached, the voltage across the test length 
of the model and across the standard resistor in series with the model 
was read on the potentiometer. The heat supplied to the model was then 
the product of the voltage drop across the test length of the model and 
the current through it. 

Model 3^ the butt -welded thermocouple, was not used in the heat- 
transfer tests. 

Independent control of the end temperatures was not possible with 
models 5 and 6; however, these models were constructed in such a manner 
that the resistance to heat flow along the models was large and the test 
length was confined to the center 3/8 inch of the models, thus minimizing 
end losses. 

The radiation heat loss from the models was determined by measuring 
the heat loss over a range of model temperatures, with no flow passing 
through the tunnel, as the test-chamber pressure was decreased toward 
zero. As it was not possible to obtain a test-chamber pressure below 
1 micron of mercury absolute, the actual radiation heat loss was obtained 
by extrapolating the measured heat loss to zero pressure. Models 1 and 
6 were later placed in a chamber in which the pressure could be main- 
tained at 0.1 micron of merciiry absolute and the results checked the 
extrapolated data. 

The room-temperature resistance of the models was checked from day 
to day and, if a resistance changed enough so that the temperature of 
the model was in doubt by more than 2^ F, the model was either recali- 
brated or replaced. 

The tests were conducted at three nominal Mach numbers, 2.0, 2.5^ 
and 3*15 and over a range of pressures from 80 to 200 microns of mercury 
absolute . 

The test gas used was comm.ercially pure dry nitrogen. 

A summary of the test data is shown in table I. 
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RESULTS AND DISCUSSION 
Temperatiire Recovery for an Insulated Body 


The usual starting point for a set of heat-transfer experiments at 
high airspeeds is the determination of the equilibrium temperature of 
the body; that is, the temperat\ire assiimed by the body in the absence of 
heat transfer. In the case of a flat plate in continuum flow, several 
analyses are available for each of the two cases of laminar and turbu- 
lent boundary layers. In these analyses, the equilibrium ten 5 >erature is 
expressed in terms of "recovery factor," r, which can be implicitly 
defined as 

Te = T^l + 2^ r 

A summary of the available experimental and theoretical results is 
given in reference l8. In brief ^ the theoretical results indicate the 
recovery factor for laminar boundary layers to be independent of both 
Mach and Reynolds numbers and to have a value of Pr^/^. All the 
turbulent -boundary-layer analyses assumed incompressible flow, hence 
Mach number effects do not appear. The analyses, in general, indicate 
that the recovery factor for turbulent boundary layers on a flat plate 
has a value close to Pr^/^ with some analyses showing a minor depend- 
ence of recovery factor on the Reynolds number. 

For the case of a cylinder with axis normal to the stream, little 
theory exists and recourse to experiment is required. Several experi- 
ments have contributed information concerning the equilibrium tempera- 
ture of transverse cylinders in high-speed air streams, notably those 
described in references 1, 19^ and 20. There is considerable disagree- 
ment among the data cited and apparently a definitive experiment with a 
systematic variation of Reynolds and Mach numbers has yet to be made. 

The equilibrium temperat\ire for the case of cylinders in free- 
molecule flow has been determined both theoretically and experimentally 
in reference 8. Good agreement was obtained between the experimental 
results and the somewhat surprising theoretical prediction that the 
cylinder equilibrium temperature would exceed the total or stagnation 
temperature of the gas stream. This phenomenon, of course, is in direct 
contrast to the corresponding phenomenon which occurs under continuum- 
flow conditions where an insulated body can have, at most, a temperature 
equal to the stream total temperature and normally does not even attain 
this temperature due to outward heat flow in the boundary layer. The 
anomaly can be explained, however, by a consideration, in the case of 
free -molecule flow, of the magnitudes of the incident and re -emitted 
molecular energy. The incident molecular energy is computed using the 
total velocity resulting from combination of the stream mass velocity 
and the random thermal velocities. When the total velocity term is 
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squared there results a term^ which is twice the scalar product of the 
vector mass velocity and the vector thermal velocity. This scalar 
product term effectively increases the apparent internal energy of the 
gas from the continuum value of 3/2 KT per molecule to a value varying 
from 2hT to 5/2 hT depending upon the speed and orientation of the 
body. The apparent internal energy becomes 5/2 IdP for high speeds and 
body angles of attack greater than zero. This is just equal to the 
internal plus potential energy per molecule of a cube of continuum gas . 
Therefore^ the energy incident on the body in free-molecule flow becomes 
equal to that of a continuum for large values of the molecular speed 
ratio for surfaces inclined at angles of attack greater than zero. The 
molecular energy which is re-emitted from the surface is assiamed to be 
equal to the energy of a stream issuing effusively from a Maxwellian 
gas, in equilibrium at some yet unspecified temperature, into a perfect 
vacuum. For the case of an insulated body, the temperature of this gas 
is that of the body. The energy of the re-emitted stream calciilated in 
this manner is equal to 2kTe per molecule. The corresponding energy 
for a continuum gas at the same temperature is 5/2 kTe per molecule if 
both thermal and potential energies are considered. For a given re- 
emitted stream temperature, a smaller amount of energy per molecule is 
transported from a body for the case of free-molecule flow than is trans- 
ported in the case of continuum flow. It is clear that, if the same total 
amount of energy, namely the incident energy, is to be removed in each 
case, the effusive stream temperature and, hence, the body temperature 
for an insulated body must be higher for free-molecule flow than for 
continuum flow. 

The equilibrium temperature for a cylinder in a free-molecule flow 
field has been shown (equation (6)) to be a function only of the mole- 
cular speed ratio, s, or the Mach number, M. (See equation (3).) 

Because of the different parameters upon which the recovery factor 
or equilibrium temperature depends in the several regimes, it is diffi- 
cult to plot all the equilibrium-temperature or recovery-factor data on 
a single graph. However, for the case of free-moleciile flow, the theory 
predicts that the equilibrium-ten^Derature ratio is a function of the 
molecular- speed ratio only. Hence, if we should plot the ratio of the 
measured equilibrium-teniperature ratio Te/T to the theoretical 
equilibrium-temperature ratio f(s)/g(s) as a function of Knudsen 
number, the deviation of the ratio from a value of unity should define 
the lower limit of free-molecule flow in terms of the Knudsen number. 

These data are shown in figure 5 and indicate that the Knudsen number 
must have a value of approximately 2 in order that fully developed free- 
molecule flow exist for a transverse cylinder. 

The recovery-factor data are shown in figure 6, plotted as a 
function of Knudsen number. It can be seen from this figure that the 
data, in the range of Knudsen numbers from 0.02 to 2.0, are correlated 
* by Knudsen number alone with no systematic Mach number effects shown. 

It is further evident that, for values of Knudsen number greater than 0.2, 
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the recovery factor exceeds a value of unity, presmably due to the 
development of free -molecule -flow effects. Also, it can be seen that 
for values of Knudsen number less than 0.2, the recovery factor has a 
constant value indicating an independence from Reynolds number effects. 

Heat-Transfer-Test Results 

The presentation of heat-transfer data for these tests is, as with 
the recovery-factor data, beset with difficulty due to lack of a common 
correlating parameter. The results of the heat-transfer tests are shown 
in one form in figure 7 in which Nuiq is plotted as a fxmction of 
^®10’ subscript lo refers to the fact that the viscosity and thermal 

conductivity were evaluated at tunnel stagnation tenj)erature, while the 
density was evaluated at free-stream conditions. The use of Nuj^g 
Rej_o suggested by Kovasznay and Tormarck, who found that this 

particular choice of conditions for evaluation of the gas properties 
eliminated the need for including the Mach number as an additional 
variable. The data of Kovasznay and Tormarck (reference l) are shown as 
a dotted line which was computed from the expression 

• Nuio = 0.580 Reio^^^ - 0.795 

which is an empirical equation resulting in the best fit of their data. 
All the heat -transfer data obtained during the present tests are corre- 
lated by the broken line with an average deviation of ^6 percent. This 
line was faired through the data by the method of least squares and may 
be represented by the relation 


Nuio = 0.132 Re 10 °*'^® 


( 12 ) 


The solid line represents equation (ll) of this paper with the accomo- 
dation coefficient a having the value of 0.9* This value of a was 
chosen to give the best fit of the data obtained in the free -molecule 
range of Knudsen numbers . 

At the other extreme range of Reynolds mombers, however, it is seen 
that the data obtained in the Ames Laboratory low-density wind tunnel 
indicate Nusselt numbers substantially below those measured by 
Kovasznay and Tormarck. When these data were obtained, it was suspected 
that the discrepancy was due to blockage effects since the model 
(l/8-inch diameter) was large compared to the diameter (l/2 inch) of that 
portion of the stream suitable for testing. Consequently, some of the 
tests were repeated using the University of California low-density wind 
tunnel (see reference 21 for description) which has a considerably larger 
(approximately 3“inch diameter) testing region. The higher Reynolds 
number data obtained in the University of California tunnel again indi- 
cated Nusselt numbers below those obtained by Kovasznay and Tormarck 
although slightly higher than those obtained in the Ames tunnel. The 
fact that the present low-density data were obtained in two different 
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wind tunnels and (in the case of the Ames test) with three different 
models^ while the data of Kovasznay and Tormarck are corroborated by the 
data of Lowell (reference 20), indicates that the difference in Nusselt 
numbers is real and is not due to experimental error. ^ 

The discrepancy may be due to one or a combination of several 
effects. First, there is the probability that a difference in tunnel 
turbulence level existed between the tests of Kovasznay and Tormarck and 
Lowell, which were conducted at relatively high pressiires, and the pre- 
sent tests which were conducted at extremely low pressures. It appears 
likely that the turbulence level in a low-density nozzle would be low due 
to the large viscous effects present.^ Another possible reason for the 
discrepancy is a difference in, the structure of the shock wave ahead of 
the cylinder. It has been shown (reference 23) that the thickness of a 
shock wave is related to the density of the stream and can be described 
in terms of a number of mean-free -molecular paths. Since the density 
at which the present tests were conducted is of the order of I/ 5 OO to 
1/1000 that of the tests of references 1 and 20, it would appear possible 
that the resultant large increase in shock -wave thickness in the present 
tests would have some effect upon the Nusselt number. 

It is also probable that a correction for the effect of viscosity 
should be applied to the readings of the static pressure tube, as was 
done in the case of the impact pressure tube readings. The value of 
this correction (which would affect the calculated Mach and Reynolds 
numbers) is unknown at present. 


C(DNCLUSI0NS 


The following conclusions can be drawn from the data for transverse 
cylinders presented in this paper: 

1. Fully developed free-molecule flow occurs for Knudsen numbers 
of approximately 2 and higher. 


^It can be observed that the heat -transfer data obtained in the University 
of California tunnel with the 0.001-inch-diameter model are in good 
agreement with the data obtained in the Ames tunnel; however, it should 
be pointed out that a shift in the temperature calibration of this 
model resulted in an uncertainty of the order of ±10 percent in the 
data obtained in the University of California tunnel. 

^In passing, it is noted that McAdams (reference 22) reported large 
changes in heat -transfer rates from cylinders due to variation of the 
turbulence level of the air stream although insufficient details of the 
original experimental work wei'e given to enable any conclusions to be 
drawn that would be pertinent to the present discussion. 
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2. For Knudsen numbers greater than 0.2, the ten5)erature -recovery 
factor exceeds iinity even though free-molecule flow is not fully 
developed. 

3. For the case of high-speed flow of a rarefied gas about cylin- 
ders, the temperature -recovery factor is primarily dependent upon 
Knudsen number . 

i)-. Over the range of conditions covered in the present tests, the 
Wusselt niJmber is a function only of the Reynolds number if the viscosity 
and thermal conductivity are based on stagnation temperature and the 
density based on free-stream conditions. 

5* In free-molecule flow, the heat-transfer data are well corre- 
lated by the theory with the accommodation coefficient equal to 0.9. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., May iL, I95I. 
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table I.- TEST DATA 


Model 

No. 

Mach 

No. 

Static 
pressure 
(Ib/sq ft) 

Stagnation 

temperature 

(°F) 

Te 

(Op) 

Tc 

(Op) 

£!£ = Tc-Te 
(°F) 

Net heat trans- 
fer from model 
(watts X 10 “®) 

Model 

No. 

Mach 

No. 

Static 
pressure 
(ib/sq ft) 

Stagnation 

ten^jerature 

(°F) 

Te 

(®F) 

Tc 

(Op) 

AT = Tc-Te 
(^) 

Net heat trans- 
fer from model 
(watts X 10 "“®) 

1 

1.96 

0.280 

69.0 

152.5 

303 

150.5 

3.03 

2 

3.19 

0.330 

63.0 

100.5 

310.3 

209.8 

25.9 


2.01 

.529 

71.0 

— 

313 

^160.5 

5.95 


3.15 

.**37 

67.7 

101.3 

308.0 

206.7 

32.7 


1.97 

.377 

68.5 

— 

308 

^161.8 

4.32 


3.12 

. 5**7 

68.5 

96.0 

307.5 

211.5 

41.0 


1.95 

.331 

68.7 

148.5 

— 

— 

— 

2 

3.17 

.378 

64.3 

105.0 

— 

— 

— 


1.99 

.418 

71.3 

148.0 

— 

— 

— 


— 

— 

— 

— 

— 

— 

— 


2.46 

.378 

71.6 

— 

300 

^ 146.5 

4.83 

3 

3.11 

.274 

70.4 

60.6 

— 

— 

— 


2.46 

.378 

71.6 

— 

400 

^ 246.5 

8.10 


3.15 

.395 

67.4 

56.5 

— 

— 

— 


2.48 

.422 

71.3 

— 

301.5 

^149.4 

5.63 


3.04 

.503 

68.2 

57.8 

— 

— 

— 


2.48 

.422 

71.3 

— 

401.5 

^249.4 

9.80 


3.17 

.546 

68.5 

57.6 

— 

— 

— 


2.49 

.468 

68.7 

— 

300 

^153.3 

6.20 


2.46 

.399 

68.5 

59.2 

— 

— 

— 


2.49 

.468 

68.7 

— 

403 

® 256.3 

10.65 


2.62 

.411 

70.5 

60.7 

— 

— 

— 


2.58 

.503 

71.0 

— 

301.5 

^152.4 

6.73 


2 . 4 o 

.250 

68.0 

60.4 

— 

— 

— 


2.58 

.503 

71.0 

— 

400 

^250.9 

11.20 

3 

2.59 

.517 

68.5 

57.3 

— 

— 

— 


2.43 

.341 

72.3 

— 

302 

^ 146.8 

5.17 


— 

— 

— 

— 

— 

— 

— 


2.43 

.341 

72.3 

— 

404 

^ 248. 8 

8.56 

4 

1.95 

.346 

68.5 

66.0 

232 

166 

106.5 


2.50 

.338 

67.5 

— 

302 

® 152.3 

4.88 


1.98 

.418 

66.0 

66.0 

223.5 

157.5 

126.8 


2.50 

.338 

68.7 

— 

404 

^252. 9 

8.00 


2.11 

.564 

65.7 

65.? 

203.5 

137.8 

132.3 


2.52 

.457 

66.7 

— 

300 

® 155.5 

6.40 


— 

— 

— 

— 

— 

— 

— 


2.52 

.457 

66.7 

— 

398 

^253.5 

10.47 

5 

1.98 

.350 

69.0 

71.5 

200 

128.5 

45.3 


2.47 

.376 

65.0 

146.5 

— 

— 

— 


2.01 

.504 

71.0 

65.5 

200.3 

134.8 

65.2 


2.59 

.493 

68.7 

144.0 

— 

— 

— 


2.06 

.239 

70.3 

79.3 

197.2 

117.9 

31.8 


2.44 

.304 

66.4 

153.0 

— 

— 

— 


— 

— 

— 

— 

— 

— 

— 


2.55 

.437 

70.3 

148.5 

— 

— 

— 

6 

2.01 

.324 

66.5 

64.1 

198.4 

134.3 

85.5 


^2.25 

.180 

72.0 

173.5 

301.5 

128 

1.88 


2.01 

.267 

68.2 

68.1 

205.4 

137.3 

76.3 


'> 2.44 

.224 

65.0 

160. 5 

300 

139.5 

2.49 


S .04 

.449 

65.0 

58.4 

176.9 

118.5 

97.8 








— 

— 

— 

— 


2.04 

.285 

68.7 

69.1 

186.1 

117.0 

67.0 

2 

2.47 

.343 

70.3 

— 

301.6 

^168 

' 17.0 


2.05 

.507 

64.3 

57.2 

193.6 

136.4 

123.3 


2.60 

.503 

67.5 

— 

305 

^191.1 

27.1 


3.18 

.296 

67.0 

59.0 

197.0 

138.0 

123.7 


2.56 

.492 

71.5 

— 

302.5 

^177.5 

21.8 


,3.36 

.332 

65.0 

53.1 

184.5 

131.4 

146.5 


2.47 

.355 

74.0 

132.7 



— 


° 2.80 

.342 

70.0 

56.4 

178.0 

111.6 

136.7 


2.39 

.304 

68.9 

136.7 

— 

— 

— 


^ 2.65 

.272 

73.0 

57.5 

181.7 

124.2 

108.2 


2.49 

.388 

69.4 

129.0 

— 

— 

— 


t> 2.49 

.217 

74.0 

62.0 

176.0 

114.0 

80.7 


2.57 

.458 

69.7 

116.0 

— 

— 

— 


^2.51 

.139 

68.5 

58.0 

174.5 

116.5 

64.9 


3.22 

.275 

73.5 

128.0 

302.8 

174.8 

19.9 


b 2.39 

.087 

72.0 

71.0 

179 

108.0 

39.6 


3.18 

.315 

66.7 

120.3 

306.5 

186.2 

24.0 


b 2.53 

.208 

72.4 

59.7 

175 

115.3 

84.0 


®Te obtained from figure 5» 

^ata taken in University of California tunnel. 
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FIGURE L- AMES LOW-DENSITY WIND TUNNEL. 
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FIGURE 2.- CUTAWAY VIEW OF POROUS NOZZLE. 



FIGURE 3.- VARIATION OF MACH NUMBER WITH DISTANCE FROM CENTER OF STREAM 
WITH AND WITHOUT BOUNDARY- LAYER SUCTION. STREAM STATIC PRESSURE 
115 MICRONS OF MERCURY ABSOLUTE. 
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FIGURE 4." DIAGRAMMATIC SKETCH OF MODEL AND SUPPORTS. 
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RATIO OF MEASURED TO THEORETICAL 
EQUILIBRIUM TEMPERATURE 
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FIGURE 5.- VARIATION OF THE RATIO OF MEASURED TO THEORETICAL 
EQUILIBRIUM TEMPERATURE WITH FREE- STREAM KNUDSEN NUMBER 
FOR TRANSVERSE CYLINDERS. 
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FIGURE 6.- VARIATION OF RECOVERY FACTOR WITH KNUDSEN NUMBER. 
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